Abstract. The majority of the well-studied 5 Scuti stars 
show frequency pairs in the power spectra with frequency 
separations less than 0.06 cd^^ (0.7/iHz) as well as ampli- 
tude variability. We examine the interpretation in terms 
of separate excited stellar pulsation modes, single modes 
with variable amplitudes, and observational problems. 
The variable-phase technique, which examines the phase 
jumps near the times of minimum amplitude of an as- 
sumed single frequency, is applied to the extensive data of 
■ ' the star BI CMi, which shows some of the most extreme 
. behavior. The following results are found for the 5 features 
<^ I in the power spectrum which could be explained as sin- 
, gle modes with variable amplitudes or as double modes: 
?H ■ for three features it can be shown that these are indeed 
^ pairs of separate pulsation modes beating with each other: 
' at times of minimum amplitude the phase jumps are ob- 
j — ■ served and both the observed amplitude and phase vari- 
fSj ations can be predicted correctly by assuming two sepa- 
rate modes of nearly equal frequencies. Artifacts caused 
by observational error, insufficient frequency resolution or 
, variable amplitudes can be ruled out. A fourth pair has a 

J ■ probable origin in two excited modes, while a 5th case is 

, inconclusive due to long time scales of variability and small 
' amplitudes. The existence of close frequencies need to be 
, taken into account in planning the lengths of earth-based 
' as well as space campaigns so that sufficient frequency 
-«.^ , resolution is obtained. Possible reasons for the existence 
' of close frequencies in S Scuti stars are considered. They 
include the dense frequency spacing caused by the pres- 
O ■ ence of mixed modes, rotational splitting as well as near- 



coincidence of the frequencies of modes with different 



^ values (the so-called Small Spacing). 
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1. Introduction 

The 6 Scuti stars are variable stars situated in the Classi- 
cal Instability Strip, which pulsate with a large number of 
simultaneously excited radial and nonradial modes, which 
makes them well-suited for asteroseismological studies. 
The amplitudes of the more dominant modes in the typi- 
cal 6 Scuti star are a few millimag, which is much higher 
than found in the Sun. It is now possible for ground-based 
telescopes to detect a large number of simultaneously ex- 
cited modes with millimag amplitudes in stars other than 
the Sun. These studies require hundreds of telescope hours 
at observatories spread around the world. 

If it can be established that 5 Scuti stars contain close 
frequencies of nonradial pulsation as close as 0.05 cd~^, 
this has important asteroseismological implications. Such 
close frequencies cannot be explained by first-order rota- 
tional splitting since almost all S Scuti stars rotate faster 
than 10 kms""'^. The existence of extreme splitting asym- 
metries, which may be produced in rapidly rotating stars 
(e. g., Goupil et al. 2000), may produce a close doublet in 
the power spectrum: the observational evidence for such 
theoretically predicted behavior is still missing. Similar 
frequencies are also predicted in the asymptotic case if 
the two modes belong to different nonradial degrees with 
£ values separated by 2. In fact, the small frequency sepa- 
rations are used as an important asteroseismological tool. 
Furthermore, it is presently debated whether the large 
number of closely spaced modes of mixed p and g charac- 
ter predicted by theory are indeed found in evolved stars 
(e.g., see Breger & Pamyatnykh 2002). Close modes might 
provide some evidence for such modes. 

It therefore appears prudent to examine whether close 
frequencies of pulsation exist in S Scuti stars and if they 
do, to look at their properties in order to determine their 
origin. 

It is important to emphasize that the presence of dou- 
ble peaks in the power spectrum of a pulsator, or the emer- 
gence of a second peak after prewhitening the first peak 
at a similar frequency, do not necessarily mean that there 
exist two pulsation modes with close frequencies. A num- 
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ber of different explanations for such a behavior need to 
be tested. These range from instrumental effects, methods 
to find the frequencies of pulsation, to other effects intrin- 
sic to the star. Before one can conclude the presence of 
closely spaced double modes, the possible problems asso- 
ciated with the available data sets need to be investigated 
in detail, in particular: 

(i) The problem of frequency resolution: A crucial re- 
quirement for the discovery of close frequencies is sufficient 
frequency resolution. Loumos & Deeming (1978) derive a 
frequency resolution of 1.5/AT, where AT is the total 
span of observations. They make the point that the gen- 
erally used figure of 1 /AT is not really correct from a the- 
oretical point of view. Furthermore, data are not equidis- 
tantly spaced and longer time gaps can occur. This can 
decrease the resolution. In an extreme case, the resolution 
of two months of observing spaced one year apart corre- 
sponds to that of a single month. Poor frequency reso- 
lution, alone, does not lead to the detection of a double 
mode. Nevertheless, the combination of poor frequency 
resolution with one of the problems listed below, can cause 
severe problems in correctly interpreting the power spec- 
trum. 

(ii) The problem of incorrect prewhitening: In the anal- 
ysis to deduce the multiple excited modes present in pho- 
tometric data, prewhitening detected frequencies is a com- 
mon tool to detect further frequencies. Prewhitening with 
an incorrect frequency value leaves a spurious peak in the 
power spectrum with a frequency close to the prewhitened 
frequency, i. e., a doublet is found which does not corre- 
spond to two separate modes in the star. Consequently, 
it is very important to examine every detected frequency 
doublet whether a single frequency with a better value can 
explain the doublet. 

(iii) Artifacts of amplitude variability of a single mode: 
Amplitude variability on the time scale of years (or even 
months) is a common feature in S Scuti stars (e.g., see 
Breger 2000a, Arentoft et al., 2001). Some 5 Scuti stars 
are known to exhibit very strong amplitude variability, 
e. g., the star 4 CVn with an average annual amplitude 
variability of 12% for the different pulsation modes and 
40 % over a decade (Breger 2000a). The techniques used 
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to detect the frequencies of pulsation usually rely on con- 
stant amplitudes. The amplitude variability shows up in 
the power spectrum as multiple peaks. Strictly speaking, 
amplitude variability associated with a single frequency 
does not lead to a double peak, but to a more complex 
structure. However, with amplitude variability the exact 
value of the pulsation frequency can sometimes not be de- 
termined, so that the problem of incorrect prewhitening 
also occurs. In some cases, with amplitude variability only 
a second peak with sufficient power is noticed, while the 
other features blend into the noise present in the power 
spectrum. 

(iv) Observational problems: These range from small 
systematic time errors in the data from one of the multiple 
observing sites to different effective wavelengths for nom- 
inally similar filters at different sites. The latter is caused 
by the fact that the amplitudes of 6 Scuti stars are ex- 
tremely wavelength dependent. These problems should be 
carefully checked for by the observers during data reduc- 
tion, since at a later stage it is difScult to discover such 
effects. 

(v) Effects of aliasing: In a multiple-mode star, alias- 
ing may result in combinations of peaks with incorrect fre- 
quencies. Prewhitening such peaks, again, leaves artifacts 
in the power spectrum, which can be interpreted as double 
modes. For small-amplitude modes, we regard this effect 
as one of the most dangerous possibilities. This needs to 
be examined carefully with statistical programs specializ- 
ing in the discrimination between various multifrequency 
solutions, such as PERIOD (Breger 1990) and PERIOD98 
(Sperl 1998). 

2. 6 Scuti stars with reported close frequencies 

A number of 6 Scuti stars have been reported to have close 
frequency pairs with separations less than 0.1 c/d. In or- 
der to estimate whether this is a common occurrence, wc 
have examined the literature for all stars listed in the 6 
Scuti star catalogue (Rodriguez et al. 2000). For most of 
these stars, very little information is available, so that the 
question of close frequencies cannot be addressed. Conse- 
quently, we selected those stars which were well-studied 
so that close frequencies (if they exist) could be found. 
Due to the inhomogeneity of the data, no rigorous, con- 
sistent definition is possible. Therefore we selected stars 
with more than 200 hours of photometric observations, 
annual observing runs of 20 days or longer, as well as 
no obvious problems of analysis, measurements, and the 
distribution of the observations in time. 17 stars could 
be selected. Of these, seven stars had reported close fre- 
quencies. This underestimates the true occurrence because 
many S Scuti stars have very small amplitudes so that only 
a few modes are observed. If we further restrict the sample 
to those stars for which more than ten frequencies of pul- 
sation had been found, every one of the seven stars had at 
least one reported close frequency pair with separations of 
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Fig. 1. The pulsation frequency spectrum of BI CMi from 

4-18 cd"'^. The double peaks in the power spectrum are 
marked with a 'd', while the mode with a highly variable 
amplitude has been indicated by a 'v'. 

less than 0.06 c/d. This is surprising, since the frequencies 
typically cover a wide range between 5 and 15 c/d. Due to 
inhomogeneity of the data and the small number of well- 
studied stars, one should be careful not to ovcrinterprct 
these numbers. A safe conclusion would be that the ma- 
jority of well-studied 5 Scuti stars show close frequency 
pairs. 

Table 1 lists the seven stars with reported close fre- 
quency pairs with separations of 0.06 c/d or less. 

Let us mention two of these stars: Alvarez et al. (1998) 
reported two close frequency pairs for the star BW Cnc 
in the Praesepe clusters. The pairs at 138.7/139.1 and 
260.6/261.5 /iHz have separations of 0.4 and 0.9 /iHz, re- 
spectively. This translates to 0.035 and 0.078 cd~^, re- 
spectively. The length of the campaign of 26d suggests a 
frequency resolution of 0.058 cd"-'^ . according to the Deem- 
ing & Loumos criterion so that at least one of the reported 
doublets has sufficient frequency resolution. Another star 
is BI CMi, which was studied extensively by the Delta 
Scuti Network for a number of years so that a frequency 
separation as small as 0.005 cd~^ could still be resolved. 

Some of these studies are quite extensive. Due to the 
large amounts of data it appears inappropriate to dismiss 
all reports of close frequencies in 6 Scuti stars as due to 
observational difficulties. It is important to examine a star 
in close detail to test the different interpretations. 

We have chosen the star BI CMi for these detailed ex- 
aminations. The choice was motivated by the fact that the 
star has extensive observations available and that we are 
intimately familiar with the data, since they were obtained 
by the Delta Scuti Network located in Vienna. 

After the pioneering study of BI CMi by Mantegazza 
& Poretti (1994), from 1997 to 2000, the Delta Scuti Net- 
work obtained more than 1000 hours of photometric data 
and deduced more than 25 frequencies of pulsation (Breger 
et al. 2002). The 1998/1999 and 1999/2000 observing sea- 
sons covered 125d and 109d. respectively. This provides a 
frequency resolution of 0.003 cd~^ for the two years. One 
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Table 1. 6 Scuti stars with reported close frequencies 



star 


Frequency 


Frequency 


Separation 


Resolution 


Sufficient? 


Reference 




{ca ) 


\ca ) 


(ca ) 


^ca ; 






BI CMi 


4.783 


4.818 


0.035 


0.003 


yes 


Breger et al. (2002) 




8.641 


8.658 


0.017 


0.003 


yes 






10.429 


10.437 


0.008 


0.003 


yes 






1 9 '^97 


1 9 'i'^n 


n 09^ 


u.uuo 


yes 




XX Pyx 


27.011 


27.102 


0.091 


0.019 


yes 


Handler et al. (2000) 




38.065 


38.110 


0.045 


0.019 


yes 




HD 18878 


11.178 


11.219 


0.041 


0.025 


yes 


Mantegazza & Poretti (1993) 


4 CVn 


5.048 


5.134 


0.086 


0.028 


yes 


Breger et al. (1999) 




6.404 


6.440 


0.036 


0.028 


yes 






6.680 


6.750 


0.070 


0.028 


yes 




FG Vir 


24.200 


24.228 


0.028 


0.039 


no 


Breger et al. (1998) 


BV Cir 


11.077 


11.128 


0.051 


0.047 


yes 


Mantegazza et al. (2001) 




12.289 


12.381 


0.092 


0.047 


yes 


Kurtz (1981) 


BW One 


11.984 


12.018 


0.035 


0.052 


no 


Alvarez ct al. (1998) 




22.516 


22.594 


0.078 


0.052 


yes 


Michel et al. (1999) 



of the important and puzzling results of the study con- 
cerns the presence of five features in the power spectrum 
which can not be explained as single pulsation modes with 
constant amplitudes. Their unusual nature modes was dis- 
covered when prewhitening the frequency of each mode 
revealed an additional peak in the power spectrum with 
a frequency very close to that of the prewhitened mode 
and/or when different amplitudes were found for differ- 
ent years. We note that the detection of these features is 
statistically very significant (amplitude signal/noise ratio 
as high as 6), so that their reality will not be questioned 
here. The frequency spectrum of BI CMi is shown in Fig. 1, 
where the location of these features is also shown. 

Before the nature of these features is examined in de- 
tails, one should eliminate the possibility of their origin 
due to observational uncertainties, complex aliasing ef- 
fects from other modes or incorrect prewhitening. While 
such problems can never be eliminated completely, except 
for feature 5 (11.26 cd'^) we consider them extremely un- 
likely as an explanation for the following reason: The data 
were subdivided into various subgroups based on year and 
observatory and a large number of new multifrequcncy so- 
lutions with a variety of hypothetical scenarios computed. 
The result was always the same: two modes with close fre- 
quencies or amplitude variability. In the case of feature 5, 
we note the existence of a pulsation mode at 10.24 cd^^ 
with a similar amplitude, so that interaction through 1 
cd~^ aliasing presents a potential influence. Again, no ob- 
vious alternate explanation could be derived, but for this 
feature we would like to advise some caution. 

Two explanations remain: pairs of close frequencies 
and amplitude variability. The next section will concen- 
trate on distinguishing between these possibilities. 



3. Close frequency pairs or just amplitude 
variability? 

In principle, the choice between the hypotheses looks sim- 
ple: one needs to consider both models and calculate op- 
timum fits of the relevant parameters to the data. The 
models with the lowest residuals between fit and observa- 
tions would then be chosen. This valid method has one 
serious drawback: in a multiperiodic star, the two models 
may lead to only slightly different residuals. Furthermore, 
the solution with the lowest residuals may not be the cor- 
rect one due to complex aliasing effects (see above) and 
the different number of free parameters fitted. 

It appears prudent to adopt a more classical, analytical 
approach. The two hypotheses both lead to a variable am- 
plitude, but to very different phasing behavior: two beat- 
ing frequencies appear as a single frequency with variable 
phase. The most extreme case occurs for two frequencies 
with the same amplitudes: there is a half cycle phase jump 
at every new cycle of the beat period. Even if the ampli- 
tudes are not identical, the phase jump still exists, but 
is progressively smeared out in time as the amplitudes 
become more and more unequal. A numerical modeling 
of the amplitude ratios and phase shifts can be seen in 
Breger (1981), where the two- frequency model could be 
ruled out for the cepheid HR 7308. 

To distinguish between the two hypotheses, we will ex- 
amine the amplitude and phase variability for each mode 
of BI CMi and compare the one- and two-frequency mod- 
els in the time, amplitude and the time, phase domains. 
All computations of amplitudes and phases of the mode 
to be tested for duplicity were carried out simultaneously 
with all the other frequencies or frequency pairs known for 
BI CMi as listed in Breger et al. (2002). The simultaneous 
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Feature 


Single-frequency model 
Frequency Amplitude in y 
(millimag) 


Double-frequency model 
Frequency Amplitude in y 
(cd~^) (millimag) 


1 


10.4314 


0.6 to 6.5 


10.4289 
10.4365 


3.6 
2.7 


2 


8.6490 


0.2 to 2.5 


8.6578 
8.6405 


1.4 
1.2 


3 


4.7826 


0.3 to 2.7 


4.7826 
4.8179 


1.3 
1.0 


4 


12.3500 


0.2 to 1.3 


12.3500 
12.3268 


0.8 
0.5 



Note that the observed amplitude range for the single-frequency model does not correspond exactly to the sums and 
differences of the amplitudes of the double-frequency model since the derived parameters are based on observed data. 



calculation eliminates any bias caused by spectral leakage 
from other modes which were assumed to have constant 
annual amplitudes. Numerical simulations showed, how- 
ever, that such bias is very minor to start with. The avail- 
able data have been obtained through the Stromgren y 
and V filters. We have used both filters to reduce the noise. 
For the analysis, the v amplitudes were multiplied by an 
experimentally determined factor of 0.632 to agree with 
the y amplitudes and any small shifts between the times 
of maximum for the two filters ignored. This is a safe pro- 
cedure because in the test we are looking for phase jumps 
of up to half a cycle. Furthermore, the results are insen- 
sitive to the chosen value for the amplitude ratio because 
of the similar coverage within each for year for the two 
colors. The data, which were obtained at different obser- 
vatories with different detectors and filters, were examined 
for systematic differences (which might lead to errors in 
the power spectra). No systematic effects or problems were 
found. 

The test relies mainly on the phase shifts of the as- 
sumed single frequency. The results can be represented 
in different graphical forms based on the uncertainties of 
the calculated points and the number of measured data 
points required for each phase range considered. Conse- 
quently, we have chosen slightly different graphical forms 
for demonstrating the results for the different frequencies 
and to strengthen the conclusions. 

Table 2 lists the fom- bc;st candidates for double modes 
in BI CMi and shows the relevant parameters for both the 
two-frequency and single-frequency hypotheses. 

3.1. The mode(s) near 10.43 cdr'^ 

The pulsation mode near 10.43 cd~^ shows strong ampli- 
tude variability on the timescale of about 132d. To apply 
the phase test we have assumed the single (best) frequency 



of 10.43142 cd~^ and computed separate solutions of the 
amplitudes and phases in 10 to 15 d intervals. This could 
be done for all three observing seasons. The separate re- 
sults were then plotted together utilizing a beat frequency 
of 132d. The results are shown in Fig. 2. 

We consider the results as surprisingly unambiguous 
and find: 

— The rapid phase change is found. The fact that the 
phase is almost constant outside the times of the rapid 
switch indicates that we have chosen a correct trial pe- 
riod. The rapid phase change is incompatible with the 
single frequency /variable amplitude hypothesis and a 
beating phenomenon of two close periods with similar 
amplitudes may be involved. 

— The fast phase changes occur near the times of min- 
imum amplitude. This is another requirement of the 
two-frequency hypothesis. 

— Both the 1998/1999 and 1999/2000 observing seasons 
show exactly the same behavior. Apart from confirm- 
ing our conclusions in favor of the two-frequency hy- 
pothesis this result indicates that the amplitudes of 
the two pulsation modes are not strongly variable on 
an annual basis. 

— The 7 nights of data from 1997 spanning 15 days 
can provide only one point in Fig. 2. This point is in 
agreement with the other data and supports the two- 
frequency model. 

— We have computed the best values of the two fre- 
quencies, amplitudes and phases by a least-squares fit. 
In fact, the two- frequency model successfully predicts 
both the observed amplitudes and phases for all three 
observing seasons. 

We conclude that the model with a single frequency 
and variable amplitudes can be ruled out from the ob- 



6 



M. Brcgcr, K. M. Bischof: Close frequency pairs in Delta Scuti stars 



1120 1140 1160 1180 1200 1220 1240 1260 1280 

8 I , , , I , , I , ^ r 




1120 1140 1160 1180 1200 1220 1240 1260 1280 

Date (HJD + 245 0000 for 1998/1999) 
(HJD + 245 0395 for 1999/2000) 

Fig. 2. Amplitude and pliase shift diagrams to test 
whether the particular mode is a single frequency with 
variable amplitude or composed of two close frequencies 
with nearly equal amplitudes. Closed circles: 1998/1999 
season, open circles: 1999/2000 season, cross: 1997. The 
rapid phase change in the middle of the diagram demon- 
strates the beating phenomenon and the existence of two 
close frequencies. 



servations and that the model with two close frequencies 
provides excellent fits to the data. 

3.2. The mode(s) near 8.65 cdr'^ 

The successful analysis for the previous (double) mode 
was repeated for the 8.65 cd~^ mode. Here, the shorter 
beat cycle of ~ 56d makes the test more difficult than in 
the previous example, because the time bins for the exam- 
ination have to be shorter in length and thereby contain 
fewer points. Fig. 3 lists the important phase information 
for the two extensive 1998/1999 and 1999/2000 observ- 
ing seasons. The data were plotted on the same diagram 
with different symbols by adopting an experimentally de- 
termined best shift of 335d (6 beat cycles). 

As for the previously examined mode, there exist reg- 
ular phase jumps of ~ half a cycle. The beat cycle is short 
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Fig. 3. Phase shift diagrams to test whether the particular 
mode at 8.65 cd~^ is a single frequency with variable am- 
plitude or composed of two close frequencies with nearly 
equal amplitudes. Closed circles: 1998/1999 season, open 
circles: 1999/2000 season. The phase jumps demonstrate 
the beating phenomenon and the existence of two close 
frequencies. 



enough for each observing season to contain two of these 
phase discontinuities. These jumps occur at minimum am- 
plitude. 

We conclude that the mode near 8.65 cd~^ is actu- 
ally double with two close frequencies and that the phase 
jumps are incompatible with the single mode/ variable am- 
plitude hypothesis. As was previously also found for the 
10.43 cd~-^ pair, the amplitudes associated with the two 
modes are constant or nearly constant on an annual basis. 

3.3. The mode(s) near 4.8 cdr'^ 

The one or two modes near 4.8 cd^^ are even more difficult 
to analyze because of the short beat cycle of 28.3 d and the 
small amplitudes near 1 mmag. Consequently, it was no 
longer feasible to examine different beat cycles separately. 
The 1998/1999 and 1999/2000 data were combined. The 
short data set from 1997 was omitted in order to avoid 
potential errors accumulated from the long time base. 

To apply the phase test we have assumed the single 
(best) frequency of 4.7826 cd~^ and computed separate 
solutions of the amplitudes and phases for each for the 20 
phase bins. 

The results are shown in Fig. 4. We can see system- 
atic shifts in the phasing of the assumed single frequency. 
These shifts correspond to the shifts calculated from the 
two-frequency hypothesis. Because of the different ampli- 
tudes associated with the two modes, a sudden phase jump 
is not expected in this case. 

We conclude that in BI CMi there exist two close fre- 
quencies near 4.8 cd~^. 
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Fig. 4. Phase shift diagram to test whether the particu- 
lar mode at 4.8 cd~^ is a single frequency with variable 
amplitude or composed of two close frequencies with un- 
equal amplitudes. The phase shift of an assumed single fre- 
quency is calculated for each phase of the beat cycle. The 
large phasing uncertainty near zero beat phase is caused 
by the very small amplitude at that phase. The diagram 
demonstrates the beating phenomenon and the existence 
of two close frequencies. 

34. The m,ode(s) near 12.35 cdr'^ 

The Fourier analyses and multiple least squares solutions 
show two frequencies at 12.3499 and 12.3268 cd~^. This 
leads to a beat period of 43d. As before, we have divided 
the data into 17 bins (each with ^ 1000 data points) 
grouped according to the phasing of the beat period. How- 
ever, because of the short beat period and the low ampli- 
tudes, we were not able find an optimum single frequency 
to apply the test to. Statistically, the best result always 
turned out to be one of the frequencies of the frequency 
pair with a variable amplitude, as would occur by the 
effect of a neighboring mode or true amplitude variabil- 
ity. We have chosen the frequency with the best residu- 
als, viz., 13.50 cd^^. Fig. 5 shows that the phasing test is 
not conclusive for this case. Both the size of the residuals 
of the multiperiodic solution and the power spectra favor 
the two-frequency hypothesis over the single- frequency hy- 
pothesis. Nevertheless, a convincing proof is lacking for 
this frequency pair. 

3.5. The mode(s) near 11.26 cdr'^ 

The mode at 11.2608 cd"'^ shows strong amplitude vari- 
ability from year to year. During the 1998/99 observing 
season, covering 120d, amplitude variability including a 
minimum was found. On the other hand, no significant 
variability was found during the llOd of the 1999/2000 
observing season. These amplitudes, assuming a single fre- 
quency, are shown in Table 3, where the formal uncertain- 
ties in the amplitudes and phasing are also shown. As 
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Fig. 5. Phase shift diagram to test whether the particular 
mode at 12.35 cd~^ is a single frequency with variable am- 
plitude or composed of two close frequencies. The phase 
shift of an assumed single frequency is calculated for each 
phase of the beat cycle. The two-frequency hypothesis pro- 
vides a slightly better fit, but the scatter of the points is 
too large to convincingly distinguish between the two hy- 
potheses. 

before, the formulae for the computation of these uncer- 
tainties assume random errors. The true uncertainties may 
therefore be somewhat larger. 

The amplitude variability indicates a long beat period. 
Extensive two-frequency modelling shows a beat period 
between 500 and 600d, or Af ~ 0.002 cd"^ between the 
two (hypothetical) modes. The modelling showed that the 
two-frequency hypothesis is indeed a strong possibility, 
but that a unique, reliable solution with two close frequen- 
cies could not be found. Including the 1991 data by Man- 
tegazza & Poretti (~ 2 mmag amplitude) did not improve 
the situation because of the long observing gap between 
1991 and 1997. 

The data in Table 3 shows why the phasing test of the 
two-frequency hypothesis has to be unreliable: the test 
relies on the phase change near minimum amplitude. This 
phase change is detected, viz. 0.12 ± 0.07 cycles, but is not 
statistically significant in the light of the above discussion 
on the uncertainties. 

We conclude that due to the long beat cycle of the 
11.26 cd~^ mode, we cannot distinguish between the 
variable-amplitude and two-frequency hypotheses. 

4. The strange case of 4 CVn 

We have seen that one of the signatures of two close modes 
is the phase shift of up to half a cycle near minimum am- 
plitude. Such a phase shift has been seen before in another 
6 Scuti star: the 7.37 cd~^ mode of 4 CVn. This mode ex- 
hibits the strongest long-term amplitude variability seen 
so far. Details can be found in Breger (2000a, b: Figs. 11 
and 12). The amplitude dropped from 15 mmag in 1974 
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Table 3. Amplitude and phase variation of the 11.2608 

mode 



Observing 


HJD 


y Amplitude 


Phasing 


Season 


246 0000+ 


(mmag) 


(cycles) 


1997 


2465-2481 


2.68 ± 0.21 


0.67 ± 0.01 


1998/1999 


3138-3170 


0.64 ± 0.15 


0.61 ± 0.04 




3171-3190 


0.34 ± 0.15 


0.53 ± 0.07 




3191-3221 


0.99 ± 0.15 


0.66 ± 0.02 




3222-3257 


1.36 ± 0.15 


0.64 ± 0.02 


1999/2000 


3514-3549 


1.68 ± 0.12 


0.62 ± 0.01 




3550-3583 


1.67 ± 0.12 


0.66 ± 0.01 




3584-3624 


1.67 ± 0.12 


0.65 ± 0.01 



to 4 mmag in 1976 and 1 mmag in 1977. In subsequent 

years, the amplitude has bc!en increasing again. The amph- 
tude variability, by itself, is not unusual: the fascinating 
change is the phase shift of 0.48 ± 0.02 cycles between 
1976 and 1977, i.e., at minimum amplitude. The determi- 
nation of accurate phases became possible when the ex- 
tensive 1996/7 data of 4 CVn led to accurate frequency 
values, which also fit the earlier data well. 

Although the half cycle shift was suggestive of a double 
mode, in our previous discussions of 4 CVn we preferred 
the hypothesis of a single mode which was re-excited at 
a new epoch. The reason was that numerical simulations 
of a double mode did not lead to an improvement of the 
fit to the photometric data relative to that of the single- 
frequency variable-amplitude hypothesis. At this stage, we 
cannot present evidence proving the correctness of either 
the single-mode or double-mode hypothesis, but note that 
in the light of the new BI CMi results, the double-mode 
hypothesis should also be kept in mind for 4 CVn. The 
reason for the less definite conclusions concerning the 7.37 
cd~^ feature in 4 CVn is the relative lack of data avail- 
able for the crucial 1974-1978 time period and the long 
time gap up to the 1983/4 data. We also note that 4 CVn 
contains other close frequency pairs such as the 6.404 and 
6.440 cd~^ pair. However, if the double-mode hypothesis 
can be shown to apply to the 7.37 cd"'^ feature, then the 
best two-frequency solution gives a separation of Af = 
0.0015 cd~^. This would then be the smallest separation 
observed in 5 Scuti stars so far. 

5. Possible reasons for close frequency pairs 

There are a number of possible explanations for the close 
frequency pairs: 

(i) Mixed modes: Theoretical models predict a very 
large number of excited nonradial modes in evolved 5 Scuti 
stars due to the dual nature of these modes (p modes in the 
envelope and g modes in the core). This dense spectrum 
of modes leads to small frequency spacing. The pulsation 
models of 4 CVn (see Breger & Pamyatnykh 2002) predict 
554 unstable modes with radial orders = to 2 over a 7 
cd~^ range. Since BI CMi is in a similar stage of evolution. 



the number of modes can be regarded as an estimate for 
this star as well. The average spacing of ~ 0.01 cd~^ is of 
the same order as the separation of the close frequencies 
foimd in this paper. If the theoretically predicted modes 
are indeed unstable, the close frequencies found may find 
a natural explanation. 

We note here that the large number of closely spaced 
modes is predicted by more than one group using indepen- 
dent codes, e. g., the models for evolved 5 Scuti stars by 
Guzik, Bradley & Templeton (2000) show a similar large 
number of predicted mixed modes. 

The problem with this explanation lies in the fact 
that the large number of theoretically predicted modes 
of mixed character have not (yet) been discovered. The 
5 Scuti stars select only a small number of the possible 
modes to pulsate with detectable amplitudes and these 
sometimes show relatively large amplitudes. It is not obvi- 
ous why the (unknown) mode selection mechanism favors 
so many close modes (but see below). 

(ii) Trapped modes: Breger & Pamyatnykh (2002) have 
explored the possibility that the modes with low kinetic 
energy are preferentially selected since they are excited 
more easily. These are the modes trapped in the outer 
envelope and the analogs of pure acoustic modes. In this 
hypothesis, only the modes with the lowest kinetic energy 
would be selected out the many possible modes. The effect 
is very strong for ^ = 1 and weaker for ^ = 2 modes. This 
scenario has been successful in explaining the observed £ = 
1 frequency spacing in the star 4 CVn and is promising for 
the t = 2 modes. Again, the arguments could also apply 
to BI CMi because of its similar evolutionary status. 

In spite of the severely reduced number of excited low- 
degree modes in the trapped mode hypothesis, two modes 
with close frequencies are still possible. Calculations of 
the kinetic energy associated with all the possible modes 
show that two almost identical modes with neighboring 
frequencies can have the same low kinetic energy. Both 
modes might, therefore, be excited. 

The average spacing of the mixed- mode hypothesis ap- 
plies here as well, so that the observed frequency separa- 
tions can be matched. However, in contrast to the mixed- 
mode hypothesis, close frequencies would be predicted to 
occur only for trapped modes, i.e., at a relatively small 
number of predictable frequencies. We are, at present, en- 
gaged in testing this hypothesis with pulsation calcula- 
tions. 

(iii) Rotational splitting: Close frequency pairs can, in 
principle, be explained by rotational splitting with only 
two of the split modes visible. A simple calculation in- 
dependent of the details of stellar models appears to rule 
out this hypothesis. The transformation from a coordinate 
system rotating with the star to the observer's coordinate 
system leads to a frequency separation near mVl, where 
Vt is the frequency of rotation. Since almost all 5 Scuti 
stars rotate faster than 10 km s"-'^, the predicted frequency 
splitting is too large. However, this calculation may be too 
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simple because of the extreme splitting asymmetries which 
can be produced in rapidly rotating stars (e. g., Goupil et 
al. 2000). Such asymmetries produce rotational splitting 
which can be very much larger or smaller than that given 
by the mfl term. An observational test of this explanation 
and subsequent theoretical modelling requires reliable ob- 
servational identification of the the m values associated 
with each mode. These are not available. 

(iv) The Small Spacing: Two modes with only a small 
frequency difference arc expected if their £ values differ 
by 2, at least in the asymptotic case (e. g., Tassoul 1980). 
This is shown by 

i^ni ~ (n + £/2 + 1/4 + a)Aiy, (1) 
where n is the radial order and £ the degree of the pul- 
sation mode, Au the (large) spacing between successive 
orders, and a depends on the phase change at the stel- 
lar surface and is, in general, a slowly varying function of 
frequency. This is related to the so-called Small Spacing, 
which is expected to be an important astcroscismological 
tool for 6 Scuti stars (see Christensen-Dalsgaard 2000). 
The fact that many £ = and 2 as well as ^ = 1 and 3 
pairs are separated by only very small frequency differ- 
ences could provide an explanation for the observed close 
frequencies. In S Scuti stars the pulsations are of lower 
radial order, so that the asymptotic case does not really 
apply. Nevertheless, models show that if this explanation 
for the close modes applies, the observed frequencies would 
provide an important input into the pulsation models. 

An important prediction of this explanation would be 
that the observed close frequency pairs consist of compo- 
nent modes with £ values differing by 2. Because of pho- 
tometric cancellation effects across the surface, the only 
pairs that we could expect to observe photometrically 
would be (0, 2) as well as (1, 3). This can be checked 
if mode identifications are available. For the close fre- 
quency pair in BI CMi with the highest amplitudes, the 
mode identifications from phase differences (see Breger et 
al 2002) suggest £ = 2 for 10.429 cd'^ and f = (or maybe 
1) for 10.437 cd~^. The available mode identifications for 
the other modes and the observed frequency spacings be- 
tween the modes (i.e., the Au value from equation 1) are 
at this stage not uniquely determined. Until more data 
become available, the A£ = 2 explanation cannot be ex- 
cluded. 

(v) Mode coupling: Frequency shifts can also be pro- 
duced by nonlinear mode coupling by modes with the same 
or even different £ values (see Buchler, Goupil & Serre 
1995, and Goupil 2000). These shifts could lead to close 
frequency pairs. This effect has had very little observa- 
tional study and additional discussion is beyond the scope 
of this paper. 

6. Conclusion 

A literature search has shown that close frequency pairs 
with separations less than 0.06 cd~^ are common among 



well-studied S Scuti stars. One of the best candidates for 
a more detailed analysis of this phenomenon is the star 
BI CMi. The variable-phase technique, which examines 
the phase jumps near the times of minimum amplitude of 
an assumed single frequency, was applied to test whether 
theses pairs are caused by separate excited stellar pul- 
sation modes, single modes with variable amplitudes, or 
observational problems. 

It was shown that at least three features are indeed 
pairs of separate pulsation modes beating with each other: 
at times of minimum amplitude the predicted phase jumps 
were observed. Both the observed amplitude and phase 
variations were predicted correctly by assuming two sep- 
arate modes of nearly equal frequencies. Artifacts caused 
by observational error, insufficient frequency resolution or 
variable amplitudes could be ruled out. A fourth pair has 
a probable origin in two excited modes, while a 5th case 
is inconclusive due to long time scales of variability and 
small amplitudes. 

The existence of close frequencies needs to be taken 
into account in planning the lengths of earth-based as well 
as space campaigns so that sufficient frequency resolution 
is obtained. 
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